Replica field theories, Painleve transcendents, and exact correlation functions 
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Exact solvability is claimed for nonlinear replica a models derived in the context of random matrix 
theories. Contrary to other approaches reported in the literature, the framework outlined does not 
rely on traditional "replica symmetry breaking" but rests on a previously unnoticed exact relation 
between replica partition functions and Painleve transcendents. While expected to be applicable 
to matrix models of arbitrary symmetries, the method is used to treat fermionic replicas for the 
Gaussian unitary ensemble (GUE) , chiral GUE (symmetry classes A and AIII in Cartan classification) 
and Ginibre's ensemble of complex non-Hermitean random matrices. Further applications are briefly 
discussed. 
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Replica field theories are notoriously known for sub- 
tleties itLvolved in performing the replica limit, n — ► 0, 
devisedEl to carve a physical observable of interest out of 
the "annealed" average 



Z n (ei 



,) = (ndet n (e fc -?i) 



(1) 



\fe=i 



often called the replica partition function; generically, 
Imefc 7^ 0. The average (• • •) runs over an ensemble 
of stochastic Hamiltonians TL which, throughput this 
Letter, will be modeled by random matricesH of pre- 
scribed symmetries. Given (|l|), spectral properties of 
TL can be obtained from the p-point Green's function 
G(ei, ■ ■ ■ , e p ) — (rife=i tr ( e k—'H)~ 1 ) for which the replica 
limit reads 



G(ei, • • • ,ep) = lim 



1 



d" 



n— >o nP dei ■ ■ ■ de p 



,<*)• (2) 



Equation (0) assumes mutual commutativity of the fol- 
lowing operations: disorder averaging, differentiation, 
the replica limit and a thermodynamic limit, if necessary. 
Being an undoubtedly correct mathematical identity un- 
der the conditions of commutativity, the recipe (Q) may 
become a pitfalloij if applied unconsciously in the context 
of replica field theories. 

A difficulty lurks behind a field theoretic prescription 
to compute the partition function (|l|) . Sketchy, an origi- 
nal random system is substituted by its \n\ identical non- 
interacting copies, or replicas. Each copy, exemplified 
by a single determinant det(e — TL), is represented by a 
functional integral over an auxiliary field which is cither 
bosonic or fermionic by nature depending on the sign of 
n. Exponentiating a random Hamiltonian TL, such a rep- 
resentation facilitates a nonperturbative averaging over 
an ensemble of stochastic Hamiltonians in (Q) and even- 
tually results in effective field theories defined on cither 
a compacifl (n > 0) or a noncompactfl (n < 0) manifold. 

The point to make is this: Since the number \n\ of 
functional field integrals involved is a positive integer, 
|n| G N, the validity of the resulting representation of 
the replica partition function Z n is restricted to \n\ G N, 



too. Unfortunately, this is not enough to perform the 
replica limit (^|) determined by the behaviour of Z n in a 
close vicinity of n = 0. Therefore, a procedure of ana- 
lytic continuation of Z n away from \n\ positive integers 
is called for. 

Kamenev-Mezard's prescription. — To keep discussion 
concrete, let us reconsider a problerrD of evaluation of 
the one-point Green's function G(e) for random matri- 
ces with Gaussian distributed entries. In what follows, 
we assume the Hamiltonian- TL to be drawn from the 
Gaussian Unitary Ensembleo (GUEjv) specified by the 
probability density Pn{TL) oc exp(— trTL 2 ) for an N x N 
complex Hermitean matrix TL to occur; N E N. The 
Green's function G(e) is determined by the replica limit 
G(e) = lhn n ->on~ 1 dZ n (e)/de with 



Z n (e) = <dct n (e-TL)) 



■HGGUEa 



n G 



(3) 



Angular brackets (• • •) stand for_a matrix integral over 
the properly normalised measurea Pn(TL)VTL. 

In the above definition of Z n , the parameter n is al- 
lowed to be an arbitrary complex number C. Therefore, 
if applied directly to (|3|), the replica limit would resultEl 
in a correct answer for G(e). This is not so, howewer, 
if one maps the partition function Z n onto a varianttl of 
the replica a modelJ, Z n i— > Z n . Following the standard 
steps of fermionic mapping, one derivesatj 



Z n {e) = (det N (ie-Q)) 



QGGUE„ 



n G 



(4) 



Contrary to the starting point (|^), the replica parameter 
n in (0) is now restricted to positive integers by deriva- 
tion. As a result, any attempt to reconstruct the Green's 
function G(e) out of Z n through the replica limit (||) will 
inevitably face the problem of analytic continuation of 
Z n away from n G N. 

Despite numerous efforts and discussions throughout 
more than two decades, no mathematically satisfactory 
idea was brought in, even though there exists a recipaJ'EJ 
to deal with the problem. Its detailed exposition can be 
found in Ref.E3; below we only recollect the facts needed 
for further discussion. 
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(i) As a prerequisite, one attempts to unveil an ex- 
plicit dependence of Z n on the replica index n which is 
only implicit in (0). To this end, the matrix integral (^) 
is evaluated approximately through a saddle point proce- 
dure which makes sense if the dimension N of the random 
matrix TC is large enough. For not too large replica pa- 
rameter n £ N (in particular, n should not scale with N), 

this yieldsEl Z^ v \e) ~ J2p=o vo1 ( G «,p) z «,p( e )> where 
vol (G n , p ) = Uj=i\PU)/ T ( n - 3 + !)] is the volume 
of Grassmannian G n ,p = U(n)/U(n — p) x U(p), and 
^n,p(e) is a known function (its explicit form is irrele- 
vant to our discussion) . The inner index p in the above 
sum Y"l f- • •) counts a hierarchy of (nonequivalent) 
causalE 2 ] saddle points contributing the integral (^|) over 
Hermitean matrix Q; inequivalence of the saddle points 
highlights a phenomenon of "replica symmetry break- 
ing" . 

(ii) Further, one seeks a proper analytic continuation 
of Z^ sp) away from n £ N. Jhe (most successful so 
far) procedure devised in Ref.Ll suggests extending the 
summation over p to infinity, J2p=oi' ") ^ E^Lo(' ' *)j 
as the group volumes vol (G„ jP ) vanish for p > 1. 
Such a proposal suffers from two major drawbacksEx (a) 
for n (fc N the group volumes grow too fast with p for 
the sum Yj£Lo(' " to converge; (b) so extended to in- 
finity the sum over p would necessarily involve a con- 
tribution from n ~ O(N) where the summand is no 



been shown by DarbouxtJ a century ago, satisfies the 
equation!] 



longer given by vol (G v 



,(e). This questions the 



self-consistency of the method as a whole making it some- 
what deficient. Despite all these drawbacks, however, the 
approach furnishestl a correct result for the large- TV GUE 
density of states in both leading and subleading orders 
in 1/JV. 

Integrable hierarchies and exact Painleve reduction. — 
On brief reflection, one has to admit that the approxi- 
mate evaluation of Z n (e) is the key point to blame for 
inconsistencies encountered in the procedure of analytic 
continuation. For this reason, we opt a route based on 
exact and, therefore, truly nonperturbative evaluation of 
the replica partition function. As improbable or fantastic 
as it sounds, this is not an impossible task. 

Our claim of exact solvability of the replica model (^) 
and the models of the same ilk rests on two observations, 
(i) To make the first, we routinely reduce Jthe average 
over Q £ GUE„ in (||) to the n-fold integralQ 

/+oo n 
Y[d\ k e- x *(\ k -ie) N Al(\). (5) 
-°° k=i 

Here, A„(A) = Y[^ >e=] IXk - \i) = det(A fc -1 ) is the Van- 
dermonde determinant!] which makes it possible to bring 
(g) to the formH 



Z„(e) =e" e r n (e;iV), n £ N, 



(6) 



which involves the Hankel determinant r n (e;N) — 
det [d* +i Ti(e;N)] ke=0 n _ v The latter, as had first 



7~n — 1 Tyi 



n > l, 



(7) 



given the initial conditions tq = 1 and f\ = e _e Z\(e)\ 
the prime stands for d/de. The structure of (0) is even- 
tually due to the [3 — 2 symmetry of the replica field 
theory encoded into A 2 in (jfj) ; (3 is Dyson's index. Equa- 
tions (g) and (0) establish a hierarchy between nonper- 
turbative replica partition functions Z n with different 
n £ N. This is an exact alternative to the approximate 
solution Zn P ^£). Equation Q), known as the Toda ln±- 
tice equatioril3 in the theory of integrable hierarchiesc3 , 
is the first indication of exact solvability hidden in replica 
field theories. ■— ■ 

(ii) The second observation borrowed from RefO con- 
cerns the fact that, miraculously, the same Toda lattice 
equation governs the behaviour of so-called r-functions 
arising in the Hamiltonian formulation of the six Painleve 
equations (PI-PVI), which are yet another fundamen- 
tal object in the theory of nonlinear integrable systems. 
Complementary to (j^), and also luckily, the Painleve 
equations contain the hierarchy (or replica) index n as a 
parameter. For this reason, they serve as a proper start- 
ing point to build a consistent analytic continuation of 
nonperturbative replica partition functions away from n 
integers. This Painleve reduction confirms exact solvabil- 
ity of replica a models and assists performing the replica 
limit (j|). 

GUE density of states revisited. — With these observa- 
tions in hand, let us consider the one-point Green's func- 
tion for GUEjv where we have a rare luxury of examining 
the replica partition function as a function of energy e, 
replica parameter n, and matrix dimension N both before 
and after replica a model mapping. 

After replica mapping, the Painleve reduction of the 
partition function Z n (e) obeying (Q) and ^ with Zo{e) = 
1 and Z\ (e) = Hn (e) [Hn is the Hermite polynomial,.-see 
mh at n = 1] materialises in the exact representationED 



Z n {t) — Z n (O)exp^y dt<pTv(t)J, neK (8) 

It involves the Painleve transcendent fw(t) — (p n (N;t) 
satisfying the Eainleve IV equation in the Jimbo-Miwa- 
Okamoto forrr£3 

K v ) 2 - 4(tv4v - <^iv) 2 + Mv&'w - 2»)(*4v + 2N ) = °- 

(9) 

The boundary condition is y>rv(*) ~ (nN/t) (l + O^ 1 )) 
as t — ► +oo. Note that (||), and therefore (||), contains 
the replica index n as a parameter. 

By derivation, Eq. (||) holds for n positive integers only 
and, generically, there is no a priori reason to expect it to 
stay valid away from n £ N. We claim, however, that it 
is legitimate to extend (||) and @, as they stand, beyond 
n £ N and consider this extension as a sought analytic 
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continuation. To prove this, we examine the partition 
function Z n (e) prior to a model mapping as given by 
(||). In the eigenvalue representation 2 !, Eq. (||) reduces 
to the TV- fold integral akin to (0), 



This is a Fredholm determinan associated with a gap 
formation probability^ within the interval (2s, +00) in 
the spectrum of an auxiliary n x n Laauerre unitacy. en- 
semble. Utilising the results of Refs.EJO we deriveEfl 



Z n (e)= / J|dA fc e- A "(A fe - e )™A^(A). 

J -°° k=\ 



(10) 



Similar to §)-(@), the appearance of A 2 N in 
( |l0| ) leads to the Hankel determinant represen- 
tation Z n (e) = Z n (e; N) = e~ Ne Tjv(e;n) with 
T N {e;n) = det [d* +l n(e; n)] fe/=Q ,„ N _ 1 - Given the 

initial conditions to = 1 and t\ = e e Z n (e; 1), the Toda 
lattice equation tn t 1 ^ — (t' n ) 2 = %h-iTn+i, with N > 1, 
follows by the Darboux theoremlij. Since n is allowed 
to be an arbitrary complex number in (^), the Toda 
equation determines a whole set of replica partition 
functions Z n (e) for all n G C. For Z n (e;0) = 1 and 
Z n (e;l) = H n (ie) (see © at JV = 1), the Painleve 
reduction of the above Toda equation readsliH 



Z n {e) = Z n (0) exp I J dt^ w (t) J , n G C. (11) 

Here, ^rv(i) = ipN(n;t) satisfies the Painleve IV equa- 
tion 

« v ) 2 - 4(t^v - Viv) 2 + 4^ v (V'w ~ 2iV)(^ v + 2n) = 

(12) 

^iv(0 ~ 



behavior 
-> +00. 

Ill) reduces to (g) be- 
ipN(n; t) that can 



and matches the asymptotic 
(nN/t) (l + 0{t- r )) at infinity i 

A brief inspection reveals that 
cause of the duality iip n (N;it)\ 
easily be verified from and ([12]) and from the bound- 
ary conditions at infinity. As the above duality between 
@ and (12) formally holds beyond n G N, we con- 
clude that (g) and ([|) considered, as they stand, at arbi- 
trary complex n furnish-.a proper analytic continuation. 
Then, it can be shownta that, in the large- N limit, the 
replica projection G(e) = lim„^o n~ 1 dZ n (e)/de of the 
so-continued partition function (^) results in the correct 
expression G(e) = e — v? — 2N for the one-point, Green's 
function G(e). The famous Wigner's semicircles for the 
level density v N (e) = -7r" 1 ImG(e) = -k~ 1 ^2N - e 2 
readily follows. 

Wigner-Dyson correlations in GUE. — The two-level 
correlation function in GUE (symmetry class A in Cartan 
classification) can be treated along the same lines. De- 
fined in terms of the two-point Green's function G(s = 
ei - ea) = G( ei ,e 2 ) as R(s) = (1/2) [Re G(s) - 1], 
it can be obtained from the replica limit G(s) = 
— linin^o n~~ 2 d 2 Z n (s)/ds 2 ; Z n (s) is the fermionic par- 
tition function. Taken at imaginary argument, it is given 
by the Verbaarschot-Zirnbauer integrals 



s J ° k=l 



Z n (is) 



(A). 



(13) 



ZJis) 



exp 



i) 



2s dt Mt)-n> + nt/2y (N| 



Here, <ry(t) = a(n;t) satisnes-|the Jimbo-Miwa-Okamoto 
form of Painleve V equationn3 

(t<j v ) 2 - (cry - ta Y )[av - t<j Y + 4ct v (ct v + n)] = (15) 



ft as t 



with the boundary conditionliil 2 ! cry (t) 
+00. 

The replica limit is governed by the behaviour of Z n (s) 
in the vicinity of n = 0. In concert with the above discus- 
sion, we assume that ( |l4| ) and (JT^) determine the desired 
analytic continuation- Expanding the solution to jl~5| ) 
around n — yieldstj ay(t) = n 2 E 2 (t) -I- 0(n 3 ), where 
E%{t) is the exponential integral En(z) = j l dt e~ zt t~ e , 
Rez > 0. As a resulty-the small-n expansion of the par- 
tition function readsll3 



\nZ n (is) 



ns 



[1 + 7 + ln(2s) 



+ E 1 (2s)-E 2 (2s)}+0(n 3 ). (16) 



Here, 7 = 0.577 ... is the Euler constant. 

To the best of our knowledge, this is the first non- 
perturbative evaluation of Z n . Implementing the replica 
limit, one derives the two-point Green's function of the 
form G(s) = 1 + 2is~ 2 sin(s)e zs . This, in turn, repro- 
duces the celebrated Wigner-Dyson two-point correlation 
functiona R(s) — 7r5(s) — s~ 2 sin 2 (s). Note that this re- 
sult holds for arbitrary s down to zero (compare to Ref.0) . 

Further examples. — The same strategy can be applied 
to demonstrate exact solvability of fermionic replica a 
models for other random matrix ensembles associated 
with the Toda lattice hierarchy. Our partial list includes 
chiral GUEEa (j&hGUE, symmetry class AIII) and Gini- 
bre's ensembleE3 of complex matrices with no further 
symmetries. A detailed account of the Painleve reduc- 
tion for these ensembles will be presented elsewherdiil 
Here we announce only small-n expansions of the cor- 
responding nonperturbative replica partition functions. 
Adopting the natation of Ref.EEl for chGUE [their Eq. 

(12) ] and of Ref.EZi for Ginibre'Sr«atrix model [their Eqs. 

(13) and (32)], we have derived^ for the chGUE 



In Z { J l \e) = n(i>bxe + J dtt[K u {t) I v (t) 

+ #„_!(«)/„+!(«)]) +0(n 2 ) (17) 

by the Painleve III reduction (l v and K v are modified. 
Bessel functions) whilst for Ginibre's complex matricesca 
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hi Z n (z, z) — n(zz) 



1 



(zz) 



N 



■ 2 F 2 {N + 1, N;N + 2, N+ 2; -zz) 



(iV + l)!(JV + l) 

lnZ ra (0,2;w,w) = n [2zz + — j - n [1 +7 + ln(uu>) + Ei(wQ) - E 2 (uu))] + 0(n d ) 



(18) 
(19) 



by the Painleve V reduction (2^2 is a hypergeometric 
function). Exact correlation functions for the above 
ensembles follow from (|l7|) - ([l^) upon implementing 
proper replica limits. Out of three, the small- n expansion 
( |l7| ) is of particular interest as it holds for arbitrary val- 
ues of tie topological charge v. Alternative calculational 
schemec3 based on the saddle point evaluation of U(n) 
matrix integrals, much in line withEl, could reproduce ex- 
act results for v half integers only, with exactness, being 
securecfij by the Duistermaat-Heckman theorenO. 

(i) The A, All I and Ginibre's random matrix mod- 
els exhaust our list of ensembles illustrating exact solv- 
ability of fcrmionic replica field theories. The integrable 
structure of all of them (as well as of those for the ma- 
trix models from B, C, and D Cartan symmetry classes 
not considered in the Letter) is related to the Toda lat- 
tice (0) whose appearance is traced back to a particular 
(3 = 2 symmetry of the corresponding replica partition 
functions. 

(ii) We expect the other random matrix ensembles ex- 
hibiting (3=1 and (3 = 4 symmetries (belonging to AI, 



BDI, CI and All, CI I, Dili classes in Cartan classifica- 
tion, respectively) to be exactly solvable as well. In those: 
cases, integrable hierarchies related to the Pfaff latticeE2l 
are likely to arise. 

(iii) Another application of the formalism developed 
would be getting further insight into controversies sur- 
rounding bosonic replicas which are known to be a total 
failure in the descriptionQ of spectral correlations in some 
ensembles-jwhile being quite successful in the description 
of othersE3. 

(iv) Finally, it would be desirable to figure out to what 
extent the nonperturbative Painleve reduction of replica 
partition functions reported in this Letter is helpful in 
the replica treatment!^! of recently advocated universal 
"zero-dimensional" random Hamiltonians which include 
interactionstS. 

I appreciate a clarifying correspondence with Peter J. 
Forrester. Craig A. Tracy is thanked for helpful conver- 
sations during a workshop on random matrices in SUNY 
at Stony Brook in February 2002. 
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